Prom the mass gap in 0{N) to the non-Borel-summabihty in 0(3) and 0(4) 

sigma-models 
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We give an analytical derivation of the mass gap of the 0{N) sigma models and investigate a 
large-order behavior of the weak coupling asymptotic expansion for the energy. For sufficiently large 
A'^ the series is sign-oscillating, which is expected from the large A'^ solution of the sigma model. 
However, for A'' = 3 and = 4 the series are sign positive. 
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INTRODUCTION 

The 0{N) sigma- model is often considered as a toy 
model for the quantum cromodynamics. It is asymp- 
totically free and dynamically generates a mass scale A 
(analog of Aqcd), although its classical formulation does 
not contain any dimensionful parameters. It is widely ac- 
cepted that the asymptotic states of the sigma model are 
the massive particles in the vector representation of the 
0{N) group. The mass of the particles should be pro- 
portional to the only mass scale of the theory: m = c A. 

The mass m is a physical quantity, while c and A de- 
pend on the regularization scheme. The coefficient c can- 
not be determined from the perturbation theory. Luckily, 
the 0{N) sigma- model can be studied nonperturbatively 
due to its intcgrability. The explicit expression for the 
coefficient c in the M S scheme was found in [H, Q : 



(1) 



We explain the idea of this calculation in the next section. 
In order to obtain c, one has to solve the integral equa- 
tion ([5]) in the weak-coupling regime at the leading and 
subleading order. At the leading order the solution was 
found by the application of the generalized Wiener-Hopf 
method The subleading order, up to our knowledge, 
was found only numerically, although with a high preci- 
sion, which allowed to conjecture the expression ([T]). The 
goal of this paper is to solve the integral equation ([5]) an- 
alytically at subleading and higher orders in a recursive 
manner. 

Another motivation for this work comes from the 
AdS/CFT correspondence where the appearance of the 
0(6) sigma- model raised the necessity of the solution of 
the equation ([5| and its generalization at next-to the sub- 
leading order \^ . 

Armed with the recursive procedure, we made an es- 
timation for the large-order behavior of the coefficients 
of expansion and found unexpected properties for = 3 
and N = 4 (see the discussion session). 



INTEGRAL EQUATION 

The idea of P, 0] was to consider the sigma-model in 
the presence of a magnetic field h coupled to a U{1) 
charge. When the value of h exceeds the mass gap, a 
finite density p of equally polarized particles is created. 
At large values of h the free energy of the system can be 
computed pcrturbatively due to the asymptotic freedom. 
Knowing the free energy density f[h], we can find the 
energy density e[p] through the Legendre transform: 



e[p] = min/i(/[/i] + ph). 



(2) 



It is convenient to introduce a running coupling con- 
stant a[/i] via the relation 
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The perturbative QFT predicts the following expansion 
for the energy density: 
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where a is evaluated at the scale p = p. 

The energy of the system in the large volume can be 
calculated also from the asymptotic Bethe Ansatz, which 
explicitly contains the mass m. The considered energy 
density e is recovered in the thermodynamical limit, in 
which the number of particles K and the length of the 
system L both go to infinity with fixed p = K / L. In this 
limit the Bethe Ansatz reduces to the integral equation 
for the density of the rapidity distribution 



x[0] - 
K[e] = 



J K[9 - 9']x[0']d9' = mcosh[6i], 9^ < B^; 

log So[9], (5) 
/c.c. . 
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The energy density and the density of particles are given 

by 



e = ni 



We see that e depends on p through the parameter B. 
To compare with the expansion (jj]) we have to consider 
the large p, or equivalently large B, asymptotics of the 
integral equation 

The integral equation (O can be rewritten in a nice 
form in terms of the resolvent for the function For 
this we first notice that the kernel K\9] can be repre- 
sented as 
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is a shift operator and {1+D 
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D^^+D^^ — . . .. This representation for the kernel can be 
easily derived if to notice the formal equality r[i6'/27r] ~ 

(l/6/)T^. 

The resolvent of x[(^] defined by 



R\0] = 



dd' 



, x[0'] 



(8) 



is analytic everywhere except on the support [—B, B] of 
x[d]- The residue of R[6] at infinity equals to 27r/9/m. 
The density distribution x[ff\ can be read from the dis- 
continuity of the resolvent on the interval [—B,B]: 



(9) 



Using (O, ([H]), and ^ we can rewrite the integral equa- 
tion ([5]) as 
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LEADING ORDER SOLUTION FOR THE 
ENERGY DENSITY 

In the following we will neglect the terms that give 
exponentially suppressed contribution to the value of e. 
In this approximation we have 
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In other words, e receives the main contribution from 
the vicinity of the branch points ±B. Therefore we will 
consider the double scaling limit 
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From ^ and (fTT|) we can express the energy density in 
terms of the inverse Laplace transform of the resolvent: 
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In the double scaling limit the equation (fTO)l is simplified: 

1 - L)2A ^ __ JJ-2A 

-R{^ + *o] - R[z - lO] = 



-nime^ e'^ , z <0. (14) 
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The inverse Laplace transform of p4)) is straightfor- 
ward: 
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To find the correct solution to p5)) we demand the fol- 
lowing analytical properties for R[s] at each order oi 1/B 
expansion: 

• R[s] is analytic everywhere except on the negative 
real axis, 

• R[s] has simple poles at s = — ?i/(2A) and simple 
zeroes at s = — 1/2— n, where n is a positive integer, 

• R[s] is expanded in negative powers of s at infinity. 

The presence of zeroes in the resolvent can be directly 
seen from the equation . The solution of the corre- 
spondent homogeneous equation should have in addition 
a zero at s = —1/2. The origin of the other stated ana- 
lytical properties are explained in the appendix. 

The most general solution of the equation ([T5)l which 
satisfies stated analyticity properties is the following: 
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The dependence of Q[s\ on B is not a consequence of 
(fTS)) and can be deduced from the considerations of the 
next section. We see that Q[s\ does not contribute to the 
leading order of the 1/ B expansion. This solution should 
be compared with the expression for G± [i^] of 0] . 
The leading order of e is given by m A$[l/2]e-^/(47r). 
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THE PARTICLE DENSITY AND SUBLEADING 
CORRECTIONS 



First note that if we apply the operator 
1 

will get 



ji-A^_j,A {D-^^'^ + to the equation ^ we 



(d'^-^/^R[0 + lO] + D''^+^/^R[e - iO]j = 0. (17) 

In the previous section we found the most general so- 
lution in the double scaling limit. Still we have to fix 
unknown coefficients Qn.m- For this we consider a dif- 
ferent regime. We take again B oo but now we will 
be interested in the values of the resolvent R[6] at the 
distances of the order B or larger from the branch points 
6 = ztB. In this case the shift operator can be expanded 
in the Taylor series D = 1 + inde — |7r^(9e)^ -t- . . .. Solv- 
ing the equation (jl7p perturbatively we can expand the 
resolvent in the inverse powers of B: 
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where e[k] = k mod 2. The perturbative meaning of 
the expansion (|18p is most easily seen in terms of the 
variable u = 9/B. The solution gives us the value 
for the particle density from the residue of the resolvent 
at infinity: 
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If we reexpand ([T5)) in the double scaling limit 
we should recover the solution obtained in the previous 
section. This condition uniquely fixes all the coefficients 

Cn,m,fc and Qn.jn- 

The expansion psp in the double scaling limit orga- 
nizes at each order of 1/B as a 1/z expansion. Therefore 
we should compare it with the Laplace transform of the 
small s expansion of R[s]. As an illustration, we give 
here the terms of these expansions which are relevant for 
calculation of the leading and the subleading orders of p 
and e: 

n _ Co. 0,0 , Ci,o,o + Ci,o,i log [^] 

8co,i,o - 3ci,o,o - 2ci,o,i + (8co,i,i + ci,o,i) log [^] 

8B^ ' 
1 e~-m = '4 ~ /Alog[#]+l+ (l-2A)log. 
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RESULTS AND DISCUSSIONS 

From the results of the previous sections we find the 
expressions for p and e at the leading and subleading 
orders: 
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Resolving the parametric dependence we obtain exactly 
a and terms in the expansion Q if a is defined as 
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Comparing (|22|) with Q we confirm the result ([T]). 

Comparing the solutions and (US]) one can find the 
higher order corrections to the energy density. At first 
four loops they are given by: 



X3 = ^, X4 - (24C(3)A2+ 



+8A2 - 42C(3)A - 28A -t- 21C(3) - 8) (23) 
X5 =-^(456C(3)A3-t-24A3-918C(3)A2 - 
-eOA^ + 609C(3)A - 140A - 105C(3) - 24) . 



The two-loop result xs coincides with the field theory cal- 
culations in [4|] . The fact that the energy can be expanded 
in power series over the running coupling constant (|22p is 
a nontrivial property of the integral equation ([5]) . This is 
a strong check of the validity of the bootstrap approach. 
This hidden renorm-group dynamics of the integral equa- 
tion was explained in 0]. 

We have found analytical expressions for Xn up to n = 
26 [sl. This allowed us to estimate the large n behavior 

of Xn- 
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For A = we have a„[0] = (-1)""^ This result is 
consistent with the fact that in the large TV limit the 
leading IR renormalon pole in the Borel plane is absent 
The large-order behavior of the coefficients Xn is 
therefore governed by the leading UV renormalon pole 
leading to the Borel-summable oscillating behavior ([24]) . 

Surprisingly. a„ ~ 1.09 for A = 1 and Un+i — 
71,^^(2.09 — 0.43(rt mod 2)), i.e. the series is non-Borel- 
summable. The Borel ambiguity is of the order A^/p^ as 
it should be from the field theory point of view (see ([4])). 

For arbitrary A > 1 the behavior of the coefficients 
a„ interpolates between those for A = 0, 1/2 and 1. 
We estimate the sign oscillation of the coefficients for 
sufficiently large n and A < A^ ~ 0.4. For 1 < A < 
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Ac all the a„ are positive. For A > 1 the asymptotic 
behavior of a„ is given by a„ ~ — airi"^A"~^, where ai 
is positive. 

The non-Borel-summability in the 0{N) sigma models 
for TV = 3 and — 4 is a quite unexpected property 
and should be understood better. For = 3 one might 
try to explain it by the presence of the instantons. For 
iV = 4 our observation is supported by the non-Borel- 
summability in the SU(N) principal chiral field model 
(PCF) at large N Q (since the 0(4) sigma-model can be 
viewed as an SU{2) PCF model). 

The method used in this paper was developed in the 
papers ^^'^ with no doubts can be applied to similar 
systems such as CP", PCF or Gross-Neveu sigma- models. 



Since T[z] is analytic everywhere except on the ray z > 
the function T[s] is analytic for Re[s] < 0. We conclude 
that R[s] is analytic everywhere in the considered region 
and has poles on the ray s < 0. 

Similarly we consider the region Re[s] < 0,lm[s] > 0. 
This explains the pole structure of the resolvent. 

Since K[9 — 9'] is analytic a,t 9 = B, is analytic 
as well. This means that expansion of R[z] at z = is 
given by R[z] = log[z] (/3o + z/3i + + ...). The inverse 
Laplace transform of this expansion gives us the stated 
behavior of R[s] at infinity. 



APPENDIX: ANALYTICAL PROPERTIES OF 
THE RESOLVENTS 

There is yet another representation of the integral 
equation ([5]). To derive it we have to extend this equa- 
tion to the whole real axis. For this we introduce a new 
function Xh[(^] (density of holes) which has the support 
complementary to [—B, B]. The hmction x/i[^] is defined 
to be such that the equation 



Note that the stated analytical properties of R[s] 
are valid only in the double scaling limit. The inverse 
Laplace transform of R[9] at finite B is an entire function 
in the s-plane. This phenomena is illustrated by the 
function f[s] = ^ — ^e~B, This is an entire function. 
But since we are interested in /[1/2] in the large B limit 
the exponentially suppressed term can be neglected and 
effectively we get a pole. 
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^mcosh[9']Q[B^ - {9'y] (25) 

is valid for any real 9'. 

Integrating the equation (|25|1 with the Cauchy kernel 
e3F' get 
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Rh[9] - T[9], Im[9] ^ 0, (26) 
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where Rh is the resolvent for Xh and T is the resolvent 
for the r.h.s of (|25)) . 

The function R[s] is defined via the inverse Laplace 
transform (see ([T3| ) and is analytic for Re[s] > 0. We 
define the analytical continuation of R[s] to Re[s] < via 
the path that does not cross the ray s < 0. To study the 
analytical properties of R[s] in the region Rc[s], Im[s] < 
we consider the equation (pS)) for Im[6'] in the double 
scaling limit and apply the integral J.^J^'^q ° ^fee^^. The 
result can be written as 
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Im[s] < 0, 
(27) 
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